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Abstract 

We check quantitatively the vahdity of some popular phenomenological approaches 
of QCD in simple models. Dispersion sum rules are considered within the ladder 
approximation of a field-theoretic model with OPE given by ordinary loop diagrams 
which are computable for any number of loops. A correlator of two currents within 
the model complies with all requirements of standard QCD sum rules approach for 
fitting low-lying resonances such as fast convergence and good stability while the exact 
spectrum contains no resonance. Optimized PT as it is inspired by the principle of 
minimal sensitivity is analyzed within a quantum mechanical model and is shown to 
work well as compared to pure asymptotic expansion in the coupling constant or Pade 
approximation. Renormalon technique is tested within another quantum mechanical 
model and is found to fail to detect its low-energy structure. 
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1 Introduction 



With considerable and steady improvement of experimental data (as a recent review see, 
e.g. requirements to the accuracy of theoretical predictions are becoming stricter and 
one encounters an urgent necessity of quantitative check of phenomenological methods used 
to treat the data {e.g. 0). In the electroweak sector of the Standard Model (SM) pertur- 
bation theory (PT) works well because of smallness of the coupling constant and there is no 
immediate need in its revision that was confirmed by the accurate estimate of the top quark 
mass based on electroweak radiative corrections 0. For processes involving strong interac- 
tions PT calculations are now of the comparable precision with data that caused efforts to 
modify ordinary PT . Also PT is known not to be complete due to existence of instantons 
and confinement. Several widely used approaches were suggested in QCD for describing 
the properties of hadrons and were successfully used for phenomenological predictions while 
their real precision has never been checked as well as their real connection with fundamen- 
tal lagrangian and properties of QCD. At present it seems that this latter problem will be 
solved only within the lattice approximation that is going to attain a status of ultimate 
judge of validity of computational methods. Not everything is yet computable within lattice 
so some phenomenological methods have no quantitative tests. In these circumstances it is 
instructive to analyze some models. 

In this paper we give a quantitative analysis of several methods used in QCD within 
exactly solvable models. They are dispersion sum rules, optimized PT, and partial resum- 
mation of PT for pinning down the non-PT effects (renormalons) . The common feature of 
the above techniques is an attempt to obtain the non-PT information on the strong interac- 
tions studying asymptotic series in strong coupling constant (within ordinary PT) or large 
momentum asymptotic expansions (within operator product expansion (OPE)). 

First we consider a two-point correlator in a field-theoretic model and restrict ourselves 
to a kind of massless (j)^ ladder approximation. The spectral density of the correlator within 
the approximation is known explicitly and does not contain any resonances. Meanwhile 
making use of the standard sum rules technique with a simple "resonance + continuum" 
model of the spectrum predicts parameters of the "resonance" very accurately in a sense 
that all necessary criteria of stability are perfectly satisfied. Though the use of sum rules 
implies that hadron properties are mainly determined by several leading terms of asymptotic 
expansion of the correlator of relevant interpolating currents in deep Euclidean domain that 
cannot be guaranteed by itself. This problem was studied within several model of quantum 
mechanics or field theories in space-time dimension less than four and sum rules proved to 
be successful though it was stressed that there is no criterion for establishing the validity 
of the technique. In two-dimensional electrodynamics (the Schwinger model) Q sum rules 
calculations were directly checked and stability of the result with respect to inclusion of 
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higher order corrections was proposed as an intrinsic criterion of sum rules apphcabihty. In 
the present paper we demonstrate within another exactly soluble model an opposite example 
when truncation of the asymptotic expansion of a correlator leads to missing some main 
properties of its spectral density while the formal stability requirement for corresponding 
sum rules is completely satisfied. This situation reflects some general property of sum rules 
and can be realized in QCD as well. This is also the first four-dimensional field-theoretical 
example with known exact answer. 

Second, within a quantum mechanical model we study the problem of resummation of 
an asymptotic perturbation series via optimization of the perturbative expansion. In most 
physically interesting models of quantum field theory the conventional perturbation theory 
forms an asymptotic series in the coupling constant which can be used for calculation of 
the Green's functions only if the effective parameter of the expansion is small enough while 
for practical purposes there exist powerful methods of refining the expansion. We also note 
that there are two different approaches to such refinement. One is more mathematical and 
consists in dealing directly with an asymptotic series applying methods of resummation like 
Borel (with hypotheses of higher order behavior) or Pade technique. Alternative approach 
relies on changing the splitting of the whole interaction into exact part and perturbation. It 
is closer to Ritz variational technique. We think that such an approach is more physically 
relevant. 

In the last part of our paper we discuss nonperturbative (power-like) corrections that 
for cases that have no simple formulation in terms of OPE are now mostly based on using 
renormalons (for a concise up-to-date review see 0). Because the expansion parameter 
- a running coupling constant - is sufficiently large for moderate energies, predictions 
differ strongly depending on a way one chooses to handle a strong coupling constant in 
infrared region. We mimic infrared renormalons of QCD by resumming a Born series for 
s-wave scattering within an exactly solvable model of quantum mechanics. Numerical esti- 
mates show that traditional technique {e.g. ^) based on introduction of nonperturbative 
power corrections fails to determine the low-energy mass scale of the model analogous to 
a typical resonance mass in QCD. An alternative approach exploiting a modified running 
coupling constant of the model and nonperturbative continuation of evolution equations into 
an infrared region gives solid and accurate estimate of this scale. 

2 Dispersion sum rules 

The problem of checking this technique is not new and fairly well understood. It is known 
that the solution to the sum rules problems is not unique. Several toy models have been 
considered where dispersion sum rules work well. We here consider an example focusing 
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mainly on the statement that the formal result is nice and satisfies all requirements of the 
technique but the exact spectral density is different and there is no intrinsic criterion to tell us 
what is going wrong. This is also a almost realistic model - four dimensional field-theoretical 
model within commonly used approximation that reveals no explicitly pathological behavior. 



2.1 The model 

We consider a field-theoretical model in four dimensions in the approximation that is equiv- 
alent to ladder massless 0^ one. In the 0^ model the expansion parameter is dimensionless 
ratio of the dimensionful coupling constant and the energy of the process in question, so 
the interaction vanishes at large energies. On the other hand at low energies the parameter 
of the expansion becomes large and the model requires a nonperturbative treatment. Thus 
it would be instructive to investigate the model along the line of ordinary QCD methods. 
Not to be taken quite seriously it gives nevertheless a way to go beyond the perturbation 
theory because here the explicit expressions for diagrams in any order of loop expansions are 
known p, ^ |T^. The model is attractive also because other attempts to break the bound of 
perturbation theory tend to be in dimensions different from 4. 

We modify the usual (p^ model slightly to make it more convenient for our purpose. The 
Lagrangian of our model reads 

£ = £0 + eiV^^A + e20^A + . . . (1) 

where Co is a free kinetic term for all fields, ellipsis stands for other interactions that are 
considered to be small and neglected, and we choose 62 = — ei = e. All questions of stability 
of the model (the existence of a stable ground state, for example) remain beyond the scope 
of our toy consideration. 

We study a correlator of two composite operators j = (pip in the ladder approximation. 
The correlator has the form 

n(g)=^|(0|Tj(x)j(0)|0)e^''-dx, n(Q2) = ^ In (^^^ + An(Q2), = -q\ (2) 

where 

where C,{z) is the Riemann's ^-function. Since the coupling constant e is dimensionful the 
expansion (j^) simulates power corrections or OPE of the ordinary QCD. Note that this 
particular subset of diagrams does not lead to infrared problems within PT as in a general 
superrenormalizable theory. Setting e^/lGvr^ = 1 we get 

p(Q2) = 167r2n(Q2) = In + Ap[Q% (4) 
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with 




MQ') = i:^[-Q^ ) I . ic(2L-i) (5) 

in analogy with QCD where the scale is given and all condensates are expressed through 
Aqcd- 

What we see first is the alternating character of the series in Euclidean direction 
(Q^ > 0) due to the special choice of interaction (|l|). This was actually the reason to create 
all decorations for the simple 0^ model in which every term has the same sign for Euclidean 
q. Note that the series (|^) is convergent for > 4 as it is seen from the asymptotic behavior 
of the coefficients at large L 

' 2l\ 4^ 



L / 



(6) 



The sum has the following closed form ||10 

oo 



n=l 

and contains no resonances. 
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1 + ^77^ -1 + 



The spectral density for the function Ap{Q ) reads 



&p{s) = -t^f^e(±-s)e(s). (7) 



Unfortunately the whole spectral density of the correlator p{s) = 1 + Ap(s) is negative in 
some domains. In case 6162 > it would have a correct positive sign but the cut would 
be situated in the wrong place of the complex plane (negative semiaxis). So, any choice of 
the interaction sign leads to unphysical spectral density and the set of ladder diagrams is 
hardly representative for the exact correlator (0|Tj(a;)j(0)|0). The same situation is realized 
in QCD. If one takes seriously the leading order correlator for vector currents, for example, 
one finds that the spectral density contains an unphysical pole at = Aq^j-, due to running 
coupling constant and does not satisfy the spectrality condition being negative at s < Aq^j-,. 
We omit these delicate points and proceed as in QCD. Namely, whether the expansion (j^) 
can be described successfully with the simple formula 

pi^^t[s) = F6{s-m^)+e{s-so) (8) 

for the spectral density p{s). 

An explicit expansion for the correlator reads 

A (n^^ 6C(3) 20C(5) 70C(7) 252C(9) 
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while the "test" form of the correlator is 

We connect expressions (§) and (|l^) by means of sum rules. There are many essentially 
equivalent ways to extract physical parameters from OPE @ and representation (|10D how- 
ever the most widely used are local duality or finite energy sum rules (FESR) [|r^ and Borel 
technique. We check them in sequence. 



2.2 Finite energy and Borel sum rules 

First we use FESR 

r%V(s)ds= r%V'^^*(s)ds (11) 
Jo Jo 

where k = 0, 1,2 because the test spectral density has three parameters to be determined. 
Eqs. ([Tl| ) now become 

F = so-6C(3), Fm2 = ^s2-20C(5), Fm^ = ^sj^ - 70C(7) (12) 

that leads to the equation for determination of the duality threshold sq 

- 2C{3)sl + 20C(5)s2 - 70C(7)so + 420C(7)C(3) - 400C(5)2 = 0. 

It has solutions sq = 16.9 and s'q = 2.24, and for corresponding parameters F = 9.67 and 
= 12.6 whilst F' = —4.97 and m'^ = 3.67. We will study the first solution because 
the second one gives an unnatural relation between the "resonance mass" and the "duality 
interval" s'q < m'^. 

Now we check the Borel sum rules approach [l^. The Borel transformation p{M'^) of the 
function piQ"^) is 

(M^\ 1 6C(3) 2005) 3507)^4209) 
the continuum contribution gives 

c(M^) = exp (-^) , (14) 
and the resonance contribution reads 

r(An = ^exp(^--j. (15) 

We have plotted these functions in Fig. 1. As we see our "test" representation accurately 
simulates the asymptotic form of the correlator at > 8. Equating the functions p(M^) 
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and r(M^) + c(M^) we obtain the Borel sum rules to determine parameters of the resonance. 
The sum rules look like ordinary QCD sum rules. Fig. 2 shows the dependence of mass 
on the Borel variable with the parameters So and F given by FESR for the different 
numbers of power corrections included. For other close values of the parameters So and F 
the results are not very different from fig. 2. A stable region is reached for 8 < Af^ < 18 
where, on the one hand, higher order power corrections are small and, on the other hand, 
the continuum contribution eq. ([T^ ) does not prevail over the resonance one eq. (|l^). The 
best stability is obtained at the optimal values So = 16.6, F = 9.39 that shows a consistency 
of the Borel sum rules and FESR. Furthermore, the curves in fig. 2 show that inclusion of 
higher order power corrections does not destroy the Borel sum rules and even enlarges the 
region of stability. 

Let us emphasize that the sum rules are perfectly saturated by the artificially introduced 
resonance and show very good stability though the exact spectral density does not contained 
any resonance singularities. The reason of this phenomenon is quite transparent: using sum 
rules we neglect the high order terms of large momentum expansion which do not affect the 
rough integral characteristics of the spectral density but are essentially responsible for its 
local behavior. The extreme sensitivity of the local form of the spectral density to the high 
orders of perturbative expansion can be easily demonstrated within considered model. For 
example, substituting ^-functions for L > 2 by units in the series (|) one neglects the terms 
with n > 2 in eq. (0) and gets the modified spectral density 



A~p{s) 



'd{A 



As we see the tiny variation of the coefficients of the asymptotic expansion leads to a drastic 
change of the spectral density at low scale s < 4/9. Indeed, the modified spectral density 
has two singular points at s = 4 and 1 while the original expression has an infinite number of 
singularities at s„ = 4/n^, n = 1, 2, ... At the same time this variation does not really affect 
the sum rules result that changes slightly (less than 10%): Sq = 17.9, F = 10.7, = 13.0. 

The lesson we draw is that our toy model containing no resonances can be well fitted by 
a standard "resonance + continuum" ansatz with good stability properties. While stability 
criteria look appropriate for determining parameters of resonance in established channels 
there is no intrinsic criterion of the reliability of QCD sum rules predictions for channels 
where there is no experimental data about the spectral density as, for instance, for gluonia. 
This makes sum rules predictions questionable. 

Our model has recently been reanalyzed in |13| where some mathematical criteria for 
validity of the "resonance + continuum" ansatz have been used in addition to a simple 
stability requirement. It was shown that making use of Holder inequalities that are valid due 
to the positiveness of the spectral density puts some additional constraints on the structure 
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of the series in and can help in ruhng out some cases of spurious resonances. Note 

however that such a technique is not apphcable in cases of non-diagonal correlators where 
the spectral density must not be positively defined. 



3 Optimized perturbation theory 

Some methods have been suggested to improve convergence of conventional perturbation 
theory. We would divide them into two groups depending on either one knows coefficients of 
expansion for ra-term that is rare or one knows the structure of interaction and makes edu- 
cated guess about the best zero order approximation. First is a pure mathematical problem 
(like a direct summation with Borel-like technique) and we don't touch it concentrating on 
the second. 

A broad class of methods is represented by general approach of (5-expansion. The basic 
idea of the optimized (5-expansion is to introduce the artificial parameter S which interpolates 
between the original theory with Hamiltonian H, and another theory, with Hamiltonian 
Ho{X) (A is a set of auxiliary parameters not present in the original theory), which is soluble 
and refiects some main properties of the theory we are interested in. One defines a new 
Hamiltonian depending on 6 

Hs = HoiX) + 6iH - Ho{X)) (16) 

and any desired quantity is evaluated as a perturbation series in S, which is set equal to 
unity at the end of the calculations. So this parameter 6 simply marks relevant terms of the 
expansion developed. Convergence of the series is achieved by an optimization procedure 
|l^ e.g. by fixing the parameters A at every finite order of the expansion in 6 according 



to principle of minimal sensitivity (PMS) at the point where the result is least sensitive to 
their variation, or the principle of fastest apparent convergence (FAC) at the point where 
the next term in the series vanishes, or some other criterion. Though the above procedure is 
not rigorous it gives good numerical results in most cases. The method has been mostly ad- 



vanced in studying the anharmonic oscillator where the convergence has been rigorously 
estabhshed ||16| . 

It seems instructive to consider the simplest model which, nevertheless, retains most 
relevant features of the real problem and, we hope, can help to gain some intuition to cure 
some difficulties of PT expansion in QCD. 
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3.1 The model 



The problem we will study is the ground state in the spectrum of the stationary Schrodinger 



equation in three dimensions 



H{a)^{r) = {-A + U{a,r))tlj{r) = Eipir), r = |r|, (17) 
where we set 2m = h = 1 and 



U{a,r) 



— ^ — a + a'^r, r < 1, 
0, r > 1, 



The choice of the potential is transparent: at a = the depth of the well is equal to the 
value at which the first s-wave bound state appears. The exact solution for the ground state 
energy E{a) is given by the transcendental equation 

y=E = «2/3Bi(eo)Ai'(6) - Ai(eo)Bi'(ei) 



Bi(eo)Ai(6)-Ai(eo)Bi(ei) ' 

eo = + « + i?) , ^1 = -a-"/'' (^^ + a - a' + Ey (19) 
where Ai(z) and Bi(z) are Airy functions [Q. Using the asymptotic expansion of Airy 



functions at large negative ^o,i (small a) we obtain an asymptotic series for the ground state 
energy 

r./ N F./ N «^ f. /3 2 \ /21 11 13\ 2 
39 41 35 ^26A 



32 87r2 Qtt^ ttV 



Numerically it reads 



E{a) = -— [1 - 1.7026a + 1.2602a2 - 0.7864a=^ + •••), (21) 

and we find that the series converges fairly badly near the point a ~ 1. In fact, the series 
diverges for any positive a because the coefficients of the expansion (pOf) in high orders grow 
factorially. That reflects the presence of a singularity of the function E{a) at the origin in 
the complex a plane. The form of the singularity can be found directly from eq. ([191) : E{a) 
has a cut along the negative semiaxis and a branch point at a = 0. For sufficiently small 
\a\ it is an analytical function for |arg(a;)| < tt; therefore the series (|20| ) is Borel recoverable 
||T9| , i.e. we can restore complete information on the function E{a) from its asymptotic 
expansion. The presence of the singularity in the Green's function of eq. (|T^ reflects the 
fact that at a = the spectrum of eq. (plT]) changes qualitatively and a discrete part of the 
spectrum appears. 
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Consider now the model ( [T7| ) within the perturbation theory framework. At r < 1 the 




potential (18) consists of two parts: one is a constant and the other depends linearly on r. 
The Schrodinger equation with the constant part of the potential only has bound states for 
any positive a with the ground state energy Eq{q.) determined by the equation 

(22) 

If a is small enough one can consider the constant part of the potential to be responsible for 
creation of the bound state while the linear term is a perturbation because it is suppressed by 
an extra power of a. Then one can search for the ground state energy of eq. ( [T7| ) as a series 
in a using perturbation theory. The ground state energy of the unperturbed Hamiltonian 
(the solution of eq. (^2])) is an analytical function at the origin of the complex a plane and 
can be expanded in a convergent series for any finite positive a 

Eaia) = — - {l-{--—:]a+{ — - — -]a^ 



4 V \2 71^ J \ 16 Gtt^ tt 
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13 


3 6 


32 ~ 


87r2 


27r4 7r6 



(23) 

However, after inclusion of perturbations the ground state energy of the whole Hamiltonian 
(the solution of eq. (|T^) ) gains the singularity and the full series (|^) becomes divergent. 



3.2 Optimized perturbative expansion 

Our purpose now is to develop the optimized perturbation theory (OPT) for our toy model. 
Following the general idea of (5-expansion we have to choose the "unperturbed" Hamiltonian. 
The Hamiltonian 

Ho{a') = -A + Uo{a',r), (24) 

where 

Uo{a\r) = \^ ' (25) 
[0, r > 1, 

is the potential of the spherical well with adjustable depth seems to be the simplest and 
most appropriate choice. Here the set of parameters A in eq. ([iBD reduces to the single 
parameter a' characterizing the depth of the well and we consider < a' < 27r^ so that the 



trial Hamiltonian (p^) has only one bound state. Then the perturbation reads 



6{H{a) - Ho{a)) = 6V{a,a ,r) = ■{ (26) 

0, r > 1. 
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Expanding in S and setting 6 = 1 we obtain the n-th order approximation for the ground 
state energy as a series 

En{a,a') = Eo{a') + E^'^\a,a') + . . . + E^'''\a,a'). (27) 



The zero order approximation Eo{a') (the ground state energy of the Hamiltonian (p^) is 
determined by transcendental equation (^2]) with parameter a' instead of a. For corrections 
we have 

E''^\a,a') = / driljQ{a',r)V{a,a',r)ilJo{a',r), 



E^'\a,a')= / / dr\roia',r)Via,a',v)Mc^',E,v)\\ (2J 

JO Eo[a') — E J 



where the quantity ipo{c('^ r) is the wave function of the ground state of the Hamiltonian ( p^ ) 
and iPe{(^' ) E, r) is the s-wave function with energy E which belongs to the continuous part 
of the spectrum of the Hamiltonian (p4|). These wave functions can be easily obtained but 
the explicit expressions are too large to be presented here. 

For arbitrary a' explicit analytical expressions for the terms of the expansion (pTf ) are 
absent and we have to analyze this expansion numerically. However, if a' is small enough 
we can expand the right hand side of eq. ( pTj) in the series in a' 

r. / n /I 2 \ , / 5 11 3 \ ,2 



a' 



c?q' / 4 /3 2 12\ , , 
+ ^ 1 + -^- 7--? + -! «' + ---K (29) 



vr^ V4 TT^ vr^y 



E(^)(a.aO = -(<^(l + ...) + a.^(i + ^)(l + 

The (5-expansion (^) generates three types of perturbative series. 

Setting a' = a in eq. ( pTf ) without further expansion of -^0(^)5 -E*-*-' (a, a) in a we reproduce 
the standard perturbation theory of quantum mechanics with the unperturbed Hamiltonian 
ifo(Q^) and the perturbation \^(a,a,r). 

Expanding Eq{oi), E^^\a,a) in a we derive the asymptotic series (pil|). We should note 
that the standard analysis of the positronium spectrum consists exactly in this machinery, 
i.e. every term of perturbative expansion around the nonrelativistic Coulomb solution is 
expanded in the fine structure constant. 

In our model if a' — a = 0(q;^) the effective parameter of the expansion is proportional 

to 

/ drr{a\r)V{a,a',r)^{a\r) ^ a. (30) 
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So En{<y, a') after expansion in a correctly reproduces the first n terms of eq. (pOl). Thus if 
we are interested in the asymptotic expansion only the choice of the initial approximation 
does not play any role. 

If, however, we intend to go beyond the asymptotic expansion in a we should choose the 
"unperturbed" Hamiltonian in order to provide the best convergence of the expansion. The 
most natural way to optimize the expansion (|27| ) is to fix the parameter a' in n-th order 
according to PMS criterion at the value a^'^^ which is defined by the equation 

or according to FAC criterion at the value a^^'" so that 

i^^"^'^U'=a-O=0. (32) 

To solve eqs. (PTI, for arbitrary a one has to proceed numerically. However, for small 
a the roots of eqs. (|3T| , |32D can be found as a series in a. This allow us to observe some 



qualitative features of OPT. At n = 1 using eq. (|29|) we find 

'1 2 
.2 

2 

.2 ■ 7r"2 



ar^ = a(l-Q + ^)« + ...) , (33) 
a(^^ = a(l-(l- + ^)a + ..) . (34) 



As n — > oo we have the formal series 

a* ~ a(l + afa + afa'^ + ...), (35) 

where 7^ stands for PMS or FAC and the coefficients a* are completely determined by 
eqs. (^, ^). Moreover, by construction 

Eoo(a,a^^^) = i?o(«^^^(«)). (36) 

Power counting arguments show that a^^'^^ — a^^'" = 0{a"'~^^), i.e. a^^^ and a^'-' have 
identical asymptotic expansions. Thus in the limit n — 00 all corrections in the optimized 
expansion vanish in our model for both the PMS and FAC optimization prescriptions and 
we obtain 

E^{a,a*) = Eo{a*{a)). (37) 

Since the function £^0(0^*) can be expanded in a convergent series in a* (eq. (|23|) ) the 
singularity of the function E{a) at a = is absorbed by the function «*(«), i.e. the 
series (|35|) must be a divergent asymptotic expansion. However, we can search for the values 
a* which satisfy eqs. ( pl| ^) at every order numerically rather than as a series in a. 

Thus we have three kinds of perturbative expansion for the ground state energy: the 
asymptotic series, the standard perturbation theory and OPT. The asymptotic series seems 
to be the most primitive tool and we expect to obtain the best results using OPT. This 
assumption is completely confirmed by numerical analysis given in the next section. 
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3.3 Numerical evaluation 



The results of the numerical analysis are given in Tables In Table |l| we present the 
exact value E{a) and the results of the asymptotic expansion En{(y) up to the n-th order 
{n = 0, 1,2,3). In Table |^ we compare the exact value E[a), the result of the optimized 
expansion i?i(a;, af ^^'^) and the results of the standard perturbation theory up to the n-th 
order (n = 0,1), i.e. the values of Eo^a) and Ei{a,a) not expanded in a. The [i,j] Pade 
approximants E^^'^^{a) with i + j <3 are shown in Table 0. 

We give numerical estimates for three different values of the parameter a that represent 
typical cases: 

• a = 0.1, the asymptotic expansion is applicable and justified; 

• a = 0.5, the asymptotic expansion is still applicable but one has to deal with high 
order corrections to achieve a satisfactory accuracy. An improvement of the asymptotic 
perturbation theory is desirable; 

• a = 1.0, the asymptotic expansion in principle can provide only ~ 2% accuracy after 
summation of ~ 13 terms. Then the terms start to grow. The asymptotic perturbation 
theory must be reformulated. 

As we can see, a naive attempt to improve the convergence of the perturbative expansion 
by simply keeping the exact (not expanded) value of En{a, a) which sums up some next- 
to-leading corrections gives a good result but is essentially insufficient if a is large enough. 
This shows that optimization is important for convergence. Though we cannot directly 
demonstrate that without optimization (for any fixed a') the series ( pTl) becomes divergent. 



as is the case for the anharmonic oscillator |16l], the numerical analysis clearly shows the 



advantage of the optimized expansion. Indeed, in all cases the best convergence is achieved 
within OPT. Even for a = 1.0 taking only the first order correction we reach ~ 2% accuracy. 
Using the PMS prescription we can also correctly estimate the error of the result. For 
example, for a = 1.0 the naive estimate is {E^^\l.O,a^^^^)/EQ{a^^'^^)y ~ 0.02 which 
coincides with the real uncertainty (see Table IH). 

We should note that in our numerical analysis for all values of a the asymptotic series (l20| ) 
is truncated far from the critical order (which depends on a) where it begins to diverge. The 
bad convergence of the series reveals itself only in the fact that its terms decrease quite slowly. 
On the other hand even the third order correction is hardly available. So the accuracy of 
the perturbation theory is restricted rather by technical reasons than by the asymptotic 
character of the series. Thus the optimization is not only a resummation prescription that is 
useful in high orders of the asymptotic expansion but also gives an opportunity to improve 
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the accuracy of perturbation theory in low orders. This is an important benefit, especially 
for nontrivial systems where high order calculations are impossible. 

A remark about the Pade approach is in order. As we can see, some Fade approximants 
are closer to the exact result than the plain asymptotic expansion. However, it is not possible 
to make a choice between various approximants until the exact result or the general structure 
of the series are known. Moreover in high orders where the asymptotic character of the series 
reveals itself the Pade theory becomes useless. The reason is that because of their specific 
structure Pade approximants can correctly reproduce only pole-like singularities while the 
function E{a) has a branch point at a = 0. On the other hand OPT seems to be the most 
appropriate tool to deal with such a problem because it speeds up the convergence of the 
perturbation series choosing the most natural initial approximation for a specific model not 
by a mathematical trick. This is an automatic summation device that does not need input 
information on the form of the singularity of the bound state energy in the coupling constant 
but reproduces it via optimization. This feature can be observed in our toy model, where 
the forms of the singularity of the ground state energy at a = in the auxiliary and original 
theories are essentially different, but OPT reproduces the correct singular a dependence 
through the optimized value of auxiliary parameter a'. This general property of optimized 
perturbation theory allows one to cope even with Borel nonsummable series []TB[ while the 
class of problems where the Pade theory can be successfully applied is quite restricted. 

Note that perturbation theory for Green's functions is much more involved and one of 
the main problems is uniformity of convergence with respect to momentum. Next section is 
devoted partly to this problem. 



4 Testing renormalon technique 



In this section we check ideas of method based on notion of infrared renormalon within the 
quantum mechanical model. The technique we further refer to as a standard one presently 
consists in resumming bubble chains with principal value prescription for singularities in the 



Borel plane. Rich phenomenology can be built on such a base pO[ though the real sensitivity 
of the approach to the infrared physics is unclear as well as an unambiguous disentangle of 



perturbative and nonperturbative (condensate) contributions I^Tl. Some other approaches 
use mostly the modified running of the coupling constant |2^, |23|, a specific recipe for scale 
setting [El, E5|, EG] or some modification of (3 function to produce evolution of running 



coupling constant at small momenta Initially there is no preference between these 

techniques because no exact results on the behavior of PT in large orders or in IR domain are 
known. Some general properties of the quantum field theory to be respected (like analyticity) 
give no much help to distinguish between possibilities. Yet in phenomenological applications. 
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the existing methods give different numerical results lying on the edge of experimental errors 
and the final selection of a working frame will eventually be based on how well the particular 
technique fits experimental data. 

In this section we investigate two different approaches within a quantum mechanical 
model that reflects some general features of renormalons. 



4.1 The model 

We consider the problem of potential scattering with 

V{r) = V^6{r - ro) (38) 



and limit ourselves to s-wave amplitudes . Such a potential can be considered as a kind 
of confining (not completely) one. We study a value of wave function at the origin (a free 
wave function is normalized to 1). The exact solution for scattering of the plain wave with 
momentum k reads 

i){k) = i){k, r = {]) = {l + Yle'^'-^ sin(fcro)) . (39) 

where we set 2m = ^ = 1 as in the previous section. To study scattering of wave packages 
with distributed momentum we consider an integral of the form 

POO 

$(A) = / ^{k)W{k,\)dk (40) 





where W{k, A) is a normalized weight function of a package depending on a set of parameters 
A, /q°° W{k, A)d/c = 1. It is more convenient to deal with a function -F(A) 

^(A) = 1 + F(A) 

so that F vanishes if the scattering potential is switched off. Because of oscillating factors 
in eq. (^) integrals (^) are not well suitable for the PT analysis (they are "Minkowskian" 
quantities). For m = r^^ > |Vo| there are no bound states in the potential (|38|) and we 
can carry out the Wick rotation because ipW is analytic in upper semi-plane (Im A; > 0) 
that corresponds to the physical sheet in energy E k"^. In "Euclidean" region the exact 
expression for il){k) becomes 

^(g)= (^1 + ^(1-6-2"/-)) (41) 

where k = iq, q > 0. The last formula can be obtained by PT from a Born series for the 
standard Lippmann-Schwinger equation of potential scattering since we deal with a finite 
range potential. Each term of the Born series for tp{q) contains contributions of different 
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kinds: a power like or "perturbative" (modeling logarithmic terms of a QCD series) and 
exponentially suppressed or "nonperturbative" (modeling power corrections in QCD). We 
are trying to use an analogy with QCD as close as possible (even through terminology) 
though one must remember a toy character of the model. The parameter Vq like Aqcd 
determines the scale at which the perturbation theory series becomes poorly convergent 
while m like p-meson mass determines the scale of nonperturbative effects. Note that we do 
not fix the sign of our parameter Vq and will study both attractive and repulsive interaction. 

At high energies {q ^ m) the expansion parameter is a running coupling constant a{q) = 
Vo/2q (trivial asymptotic freedom as in superrenormalizable theories) 

^{q) = r'{a) + rnm,Vo,q) 

where 

oo 
n=0 

and ip^^{m, Vq, q) stands for exponentially suppressed "nonperturbative" terms 

r^(m,\/o,g) = ^e-2''/- + ... (43) 
Zq 

Within the present model we classify terms with respect to their behavior at large q: power 
like vanishing - PT, faster than any power - non-PT. With the standard renormalization 
group (RG) terminology we have 

/?(«) = <1^ = (44) 

Resummation in eq. ( ^2]) (in the spirit of RG) results in definition of a new (renormalization 
group improved) "running" coupling constant 

- ^ (45) 

with a j3 function 

P^'i^as) = = -aasil - a^s). (46) 

The use of this expansion parameter allows us to improve the perturbation theory and to 
sum up all "perturbative" power terms of the series (^2]) 



r'{q) = l + aUq). (47) 

Now we turn to consideration of the wave package of a specific form given by the following 
weight function 



W{q,Q) = Q 



This weight function has a bump of the width ^/SQ located at g ~ (5/2 so the above wave 
package can be considered as a "probe" of the scattering potential at the scale ~ 2/(5. 
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4.2 Borel resummation and modified f3 function 



It is easy to see that our observable pO|) suffers from tlie renormalon. Substituting ^jj{q) in 
eq. (HOf) by its asymptotic expansion (^Tf) we obtain 



/•oo 

F-%a)= / aUq)W{q,Q)dq. (4J 
Jo 



The quantity lias a typical structure of QCD observable containing renormalon, i.e. 

it is an integral of some weight function multiplied by a running coupling constant over the 
interval that includes strong coupling domain. Note that in our model the use of the running 
coupling constant (|^) in the integrand accounts for all perturbative corrections. One can 
only dream about that in QCD where a trick based on a pure assumption motivated by the 
"naive nonabelianization" is used. 

After integration we get the series with factorially growing coefficients 

oo 

= (49) 

n=l 

Properties of the series (^) depend crucially on the sign of a or Vq. Let us consider first 
repulsive potential Vq > 0. Then the alternating sign series (^9]) is Borel summable (in QCD 
it corresponds to an ultraviolet renormalon). The Borel image 

~ ?/ 

F^^iu) = --— (50) 
I + u 

has a pole at m = —1 and is a regular function on the positive semiaxis. So the Borel 
summation leads to an unambiguous result (a = Vo/2Q) 

-l/a / X 

FBia) = El - -1 (51) 



where Ei(x) is the integral exponent ||T8 



oo 



E,{x) = / —dt. 

Jx t 

One could naively hope that the Borel summation provides us with nonperturbative informa- 
tion on the function F{Q) and on the magnitude of nonperturbative parameter m. However 
numerical analysis shows that Fb{Q) does not approximate well the exact function F{Q) 
for intermediate Q m (Fig. 3). Indeed, aas{q) is the best expansion parameter (exponen- 
tially accurate) at large momenta but ip°-^{q) does not approximate well the function il){q) 
at small q and eq. (|5lD tells us nothing about the parameter m that measures exponentially 
suppressed "nonperturbative" contributions. Alternatively, one can compute the function 
F{Q) within the modified perturbation theory with sufficient accuracy even at very small 
Q. For this purpose one has to choose a relevant expansion parameter. First we write 

aas{<l) = a^\q) (l - ^««(g))" = ««(?) E (^)^«(?)", 
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where /i is a parameter. Now we limit ourselves to only two terms of this expansion that is 
reasonable in PT region and define new expansion parameter 

with a (3 function 

that is a PT transformation. The perturbation theory series for ^p{q) in a^^-* reads 

<^(g) = l + af (g)+0(a(gf). (53) 



Taking the first order term in eq. (|53D and fixing the parameter /i at some value we find the 
function F^{Q) = Oi^^KQ)W{q, Q)dq to be very close to the exact function F{Q) up to the 
very small Q (Fig. 3). Note that at very large Q the function Fb{Q) becomes closer to exact 
result F{Q) than F'^{Q) because Fb{Q) and F{Q) have the same asymptotic expansion. On 
the other hand lifting this too strong condition we get the function F^[Q) that approximates 
the exact function F{Q) uniformly for all finite Q. Because at the optimal the function 
ip^'^\q) is close to ip{q) for all finite q this approximation is universal in a sense that it works 
well for various forms of scattering package. The optimal value /i°^*(m, Vq) can be extracted 
from experiment (in our model the exact solution plays the role of experimental data). It 
turns out to be very sensitive to the variation of m. So in this way we obtain the real 
information on the scale m at which nonperturbative effects become important. 

The case of attractive potential V{r) = —Vo6{r — ro), Vq > at first sight seems to be 
completely different. The running coupling constant 

aas{q) 



2q-Vo 

with P function /3(a) = —a{l + a) has a "Landau pole" at g = Vq/2 > 0. The series 
now becomes 

oo 

F'^«(a) = ^n!a" (54) 

n=l 

and it is not Borel summable: we encounter an "infrared renormalon". Its Borel image has a 
pole on the positive semiaxis in full analogy with QCD (when the QCD (3 function is taken 
in one loop approximation). So Borel procedure leads to an ill-defined representation in case 
of attractive potential. Standard technique could consider it as a signal of the presence of 
nonperturbative contributions. One should stress however that exact function undergoes 
no qualitative change in the low energy domain (see also I^O])- Following the line of QCD 



18 



renormalon technique we define the result of Borel summation by deforming the integration 
contour in the complex u plane. The result obtained in this way depends on the specific 
form of an integration contour while an appropriate nonperturbative part (condensates or 
power corrections in QCD) must cancel this dependence. In our model we use the principal 
value (PV) prescription to define the sum of the series ( ^4]) 



FM = PV r e-"/"-^^ = ^Ei f i) - 1 (55) 
JO 1 — u a a \aj 



where Ei(x) is an integral exponent 



Ei{x) = PV / -dt, 

J —oo t 

= Vo/2Q > 0. Within the renormalon picture one should search for the exact function 
F[Q) in the form 

F{Q) = Feia) + Ce-^/" + . . . (56) 

Here the first term has the same asymptotic expansion as the exact function F, a constant 
C gives the leading exponentially suppressed correction and eUipsis stands for "higher twist" 
contributions. The power of the exponent in (^) is determined by the position of the pole 
in the Borel image. To find the value of C one has to use purely nonperturbative method or 
extract it from experiment. Note that the form of first term, the value of the constant C, and 
high order corrections do depend on summation prescription while the whole sum does not 
by construction. The parameter C{m, Vq) is quite sensitive to the variation of m so it can be 
considered as a "probe" of nonperturbative effects what is the main issue of the renormalon 
calculations. However, we find (Fig. 4 that eq. ( ^61) poorly approximates the exact function 
F[Q) for Q m. The reason is the same as in the case of Borel summable series. Namely, 
aas{q) is a bad expansion parameter at small momenta. The only difference with previous 
case is that for Borel nonsummable series the running coupling aas becomes singular at some 
point and in principle can not be used for approximation of regular function '4){q). Stress 
once again however that it still accumulates all PT terms exactly as in the previous case. 
Moreover eq. (|5^) has no relation to the exact function F{Q) because the leading non-PT 
term in the expansion of F{Q) at large Q has completely different a dependence 




As a consequence, the parameterization (pBD is not universal i.e. one gets essentially different 
values of C for different scattering packages. 

Again, introducing a running coupling constant of the form (|5^ with an appropriate /i 
we obtain a uniform approximation of the function F{Q) for all Q (Fig. 4). In QCD this 
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prescription corresponds to the use of (probably mass dependent) RG equation for strong 
coupling constant with infrared regular solution 

Note that in this way as within renormalon picture we determine the scale where non- 
perturbative effects become crucial rather than find the exact form of F{Q). However, the 
use of an appropriate expansion parameter allows us to achieve much higher accuracy at 
intermediate momenta in calculations of the integral observables like F{Q). 

5 Conclusion 

While it is clear that sum rules should be applied with great care to cases where the physical 
spectrum is not known we found that there is no intrinsic criterion within the technique 
for detection that something goes wrong. All criteria of convergence, stability, existence of 
a mass gap are perfectly satisfied. Still this is not completely artificial construction but a 
model based on an approximation of quantum field theory. 

The optimized perturbation theory shows that knowledge of the potential allows to re- 
duce the necessary number of PT terms to a bare minimum still providing good accuracy. 
We would like to stress that in our model we used more than simple play with a finite 
piece of PT series, we chose the leading approximation differently. In QCD, however, the 
choice of zero order approximation is practically unique because there is no other solvable 
approximation but free relativistic field theory and no real optimization is therefore possible 
at present. In this respect a simple change of a renormalization prescription has no deep 
physical justification though can be quite successful phenomeno logically in some particular 
cases. 

The analysis of two possible ways of extracting information on low energy domain of a 
quantum mechanical model shows that approach based on optimization of the PT through 
introduction of a flexible expansion parameter using the freedom of choice of the scheme 
is more efficient for moderate energies than direct Borel resummation technique. In our 
model also the singularity of Borel image for attractive potential does not correspond to 
leading non-PT asymptotics of the exact function that is one of main reasons for using the 
renormalon technique in phenomenological applications of QCD. Though obtained in a toy 
model, these observations may serve as a ground for using a modified running of the coupling 
constant of QCD in the infrared domain for phenomenological applications instead (or in 
addition to) the standard renormalon technique. 
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Tables 



a 


E{a) 


^o(a) 


Ei{a) 


E2{a) 


E^{a) 


0.1 


0.2101 


0.2500 


0.2074 


0.2105 


0.2104 


0.5 


2.426 


6.250 


0.9292 


2.898 


2.284 


1.0 


2.144 


25.00 


-17.57 


13.94 


-5.721 



Table 1: The exact ground state energy E{a) and the result of the asymptotic expansion 
En{a) up to the n-th order (n = 0, 1, 2, 3) (in units of 10~^). 



a 


E{a) 




E^{a) 


Ei{a^ a) 


0.1 


0.2101 


0.2104 


0.2430 


0.2092 


0.5 


2.426 


2.424 


5.531 


1.915 


1.0 


2.144 


2.116 


20.13 


-4.200 



Table 2: The exact ground state energy E[a), the result of the optimized expansion up to 
the first order £'i(q;, af^'^) and the results of the standard perturbation theory up to the 
n-th order En{a, a) {n — 0, 1) (in units of 10~^). 



a 


E{a) 


^[i'il(«) 


^[0.2] (a) 




E[i'2l(a) 


^[0'3l(a) 


0.1 


0.2101 


0.2104 


0.1874 


0.2104 


0.2104 


0.2104 


0.5 


2.426 


2.366 


2.340 


2.430 


2.424 


2.562 


1.0 


2.144 


0.5388 


5.759 


1.833 


1.600 


4.331 



Table 3: The exact ground state energy E{a) and various [i,j] Pade approximants E^'^'^^{a) 
with i+ j <3 {in units of lO"^). 
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Figure Captions 



Fig. 1. The resonance contribution r{M'^) to the Borel sum rules and the functioin f{M'^) ~ 
p{M^)-c{M^). 

Fig. 2. The mass plotted as a function of Borel variable with two (a), three (b) and 
four (c) orders of power corrections included. The arrows mark stability interval of the Borel 
sum rules. 

Fig. 3. Numerical results for repulsive potential Vq — 1, m — 3. 
Function F^{Q)/F{Q), the optimal value // = 6 (line a). 
Function Fb{Q)/F{Q) (line h). 

Fig. 4. Numerical results for atractive potential VQ — l,m — ?>. 
Function F^^{Q)/F{Q) for the optimal value // = 3.5 (line a). 
Function Fb{Q)/F{Q) (line h). 

Function {Fb{Q) + Ce''^/'^) / F{Q) for an optimal value C = -0.06 (line c). 
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Fig. 1 
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Fig. 4 
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